ftn725998 


CONTRACT  REPORT  NO.  37 
3SR-201 


HELP  CODE  SOLUTIONS  TO  TWO  TEST 
PROBLEMS  IN  ARMOR  PENETRATION 
TOPICAL  REPORT 


V 

I 

f 


Prepared  by 

Systems,  S^clence  and  Software 
La  Jolla,  California 


May  1971 


D  D  C 

lb  JUL  »  S"  1 
D. 


This  docuwnt  has  been  approveti  for  public  release  and  sale: 
Its  distribution  Is  unllelted. 


NATIONALTEcflNlCAL 

i.ici"kDkAATION  SERVICE 


U.S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 

BALLISTIC  RESEARCH  LABORATORIES 

ABERDEEN  PROVING  GROUND,  MARYLAND 


llc< 


Unclassified _ 

DOCUMENT  CONTROL  DATA  .  R  &  D 

(S^eurltf  etmBtHU  mlton  mi  lltl*,  b^dy  ul  ^bttrmt  t  and  twI^Mlng  mmptmtion  wmf  b»  #nr»f»rf  wh9n  >h»  owrmll  npptt  It  clm§9Ui0^j 


•  HI  PONT  Tl  TLI. 

HELP  Code  Solutions  to  TVo  Test  Problems  in  Armor  Penetration 


10.  OlOtniVUTION  tTATCMCNT 

This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


II.  4atTn*CT  : 

TVo  test  problems,  each  the  impact  of  a  right  circular  cylinder  on  a  target  plate 
of  finite  thickness,  at  a  velocity  in  the  ballistic  range,  are  solved  using  an  Eulerian 
computer  program,  HELP,  for  the  solution  of  axisymmetric  time- dependent  flows.  Test 
Problem  Number  1  is  an  iron-on-iron  impact  and  Test  Problem  Number  2  is  an  iron  pro¬ 
jectile  impacting  aluminum.  In  all  cases  the  materials  are  treated  as  compressible  and 
deviatoric  stresses  are  obtained  using  an  elastic-plastic  constitutive  equation. 

Computed  results  from  these  first  two  HELP  problems  are  given  in  considerable 
detail  in  order  that  quantitative  comparisons  of  these  solutions  can  be  made  with 
Lagrangian  calculations  being  performed  by  other  research  teams. 

\ 


ta.  tPONaoRiNa  militamv  activitv 

JSA  Aberdeen  Research  §  Development  Center 
Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground.  Maryland  21005 


II.  lU^flLEMCNTARV  NOTE! 


V 

\ 


*.  RLHONT  DATE 

May  1971 

M.  CONTRACT  OR  ORANT  NO. 

DAAG  07-68-C-0931 

*.  PROJECT  NO. 

RDT8E  1T650212D620 


-  w. 


7m.  TOTAL  NO.  Or  PAQCt 

76.  NO.  or  Ncrs 

163 

3 

M.  ORISINATOR’S  REPORT  NUnniERIII 


3SR-201 


M,  OTHER  REPORT  NOIII  Mnir  otfiRT  numb»n  Mai  may  6*  maI#ip4 
Mf«  nport) 

BRL  Contract  Report  No.  37 


U,  REPORT  fSCURITV  CLAISIPIC ATION 

Unclassified 

ab.  cMou^ 


I  OHiCiNAriNO  activity  (CbfjbOMl*  AUlflorj 

Systems,  Science  and  Software 
P.O.  Box  1620 

La  Jolla,  California  92123 


Unclassified 

^•curity  fcl»««lftc«tioii 


Impact 

Armor  Penetration 
Numerical  Solutions 
Material  Response 


Unclassified 

8*curity  Clasaincatioa 


ballistic  research  laboratories 

CONTFACT  REPORT  NO  37 
3SR-201 

MAY  1971 


HELP  CODE  SOLUTIONS  TO  TWO  TEST 
problems  IN  ARMOR  PENETRATION 

TOPICAL  REPORT 


L.  J.  Hageman 
J.  M.  Walsh 

Prepared  By: 

SYSTEMS,  SCIENCE  AND  SOFTWARE 
La  Jolla,  California 


This  document  has  been  approved  for  public  release  and  sale; 
Its  distribution  is  unlimited. 


Contract  No.  DAAG  07-68-C-0931 
and 

RDT5E  Project  No.  1T650212D620 


PROVING  GROUND, 


ABERDEEN 


MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


CONTRACT  REPORT  NO.  37 
3SR-201 


» 


UHageman/JMWalsh 
Aberdeen  Proving  Ground,  Md. 
May  1971 


HELP  CODE  SOLUTIONS  TO  TWO  TEST 
PROBLEMS  IN  ARMOR  PENETRATION 

TOPICAL  REPORT 


ABSTRACT 

Two  test  problems,  each  the  impact  of  a  right  circular  cylinder  on 
a  target  plate  of  finite  thickness,  at  a  velocity  in  the  ballistic  range, 
are  solved  using  an  Eulerian  computer  program,  HELP,  for  the  solution  of 
axisymmetric  time- dependent  flows.  Test  Problem  Number  1  is  an  iron-on- 
iron  impact  and  Test  Problem  Number  2  is  an  iron  projectile  impacting 
aluminum.  In  all  cases  the  materials  are  treated  as  compressible  and 
deviatoric  stresses  are  obtained  using  an  elastic-plastic  constitutive 
equation. 

Computed  results  from  these  first  two  HELP  problems  are  given  in 
considerable  detail  in  order  that  quantitative  con^arisons  of  these 
solutions  can  be  made  with  Lagrangian  calculations  being  performed  by 
other  research  teams. 
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I.  INTRODUCTION 


General  objectives  which  have  been  assigned  to  S®  in 
the  present  armor  penetration  program  consist  of: 

(1)  The  development  of  certain  advanced  Eulerian 
numerical  techniques  for  the  solution  of 
problems  in  armor  penetration. 

(2)  The  application  of  these  techniques  to  a 
range  of  problems.  These  applications  in¬ 
corporate  accurate  material  descriptions 
and  configurations  of  practical  concern  in 
an  effort  to  provide  a  realistic  theoretical 
description  of  tV  -rmor  penetration  process. 

(3)  The  solution  of  two  specified  test  problems. 

These  test  problems  are  not  intended  to  re¬ 
present  realistic  armor  penetrations,  since 
use  is  made  of  simple  highly- idealized 
material  properties  and  projectiles  that 
are  right  circular  cylinders.  The  objective 
of  this  work,  requested  by  the  Joint  Techni¬ 
cal  Coordinating  Group  for  Munitions  Effec¬ 
tiveness  (JTCG/ME) ,  is  to  provide  solutions 

to  problems  which  are  very  well  defined  so  that 
meaningful  comparisons  can  be  made  between 
the  Eulerian  solutions  by  S®  and  the  Lagrangian 
results  by  other  research  groups  participating 
in  the  JTCG/ME  program. 

The  purpose  of  the  present  paper  is  to  report  the 
Eulerian  solutions  to  the  problems  assigned  by  the  JTCG/ME, 
Item  (3)  above.  In  Section  II  below,  the  problems  are  stated 
and  the  requested  information  is  listed.  In  Section  III  the 
solutions  are  given.  All  of  the  information  requested  by  the 
JTCG/ME  is  presented  in  that  section. 
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Both  of  the  solutions  given  in  Section  III  were  com¬ 
puted  using  the  HELP  code,  a  hydrodynamic  elastic-plastic 
multi -material  program  which  is  being  developed  and  tested 
under  Task  (1)  above.  The  HELP  code  used  as  its  starting 
point  the  code  RPM,  a  compressible  hydrodynamic  rigid-plastic 
one-material  Eulerian  code  which  was  developed  under  auspices 
of  the  Ballistic  Research  Laboratories  and  which  is  used 
extensively  in  S*  programs. 

A  considerable  amount  of  related  research,  not  required 
to  meet  the  JTCG/ME  specifications,  will  be  mentioned  briefly 
but  will  not  be  detailed  in  the  present  report.  This  work 
consists  of  Eulerian  solutions  obtained  using  alternative 
physical  models  and/or  Eulerian  programs.  One  of  these  is  an 
incompressible  rigid-plastic  solution  to  Test  Problem  No.  1 
(a  like-material  impact).  This  solution  made  use  of  the  SLUG 
code,  a  program  developed  as  part  of  this  effort  by  adding  a 
rigid-plastic  capability  to  the  S®  incompressible  Eulerian 
code  SURGE.  The  SLUG  calculation  has  been  described  in  pre¬ 
vious  reports  and  discussed  at  JTCG/ME  meetings.  Also,  a 
solution  to  Test  Problem  No.  2  (an  unlike-material  impact) 
using  the  DORF  code,  was  obtained  and  reported  in  early  1969. 
DORF  is  a  two-material  compressible  hydrodynamic  code  with 
strength  included  in  the  rigid-plastic  approximation.  Results 
from  these  and  other  Eulerian  variations  of  the  test  problems 
are,  in  most  respects,  in  substantial  agreement  with  the  HELP 
solutions  given  in  Section  III.  Some  differences,  c.g.,  the 
absence  of  tensions  in  the  incompressible  model,  have  been 
noted  previously.  The  HELP  solutions  are  presented  in 
Section  III,  however,  because  they  include  the  full  range  of 
physical  phenomena  (compressibility,  an  elastic  as  well  as  a 
plastic  regime,  and  sharp  Eulerian  interfaces)  which  occur  in 
the  problems  as  stated  by  the  JTCG/ME. 
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II.  STATEMENT  OF  PROBLEMS 

A.  INITIAL  CONDITIONS 
1 .  Test  Problem  No.  1 

Test  Problem  No.  1  is  an  axisymmetric  iron-on-iron 
impact  in  which  a  right  circular  cylinder  (radius  0.375  cm, 
length  1.15  cm)  impinges  at  9.0  x  10“  cm/sec  on  a  flat  plate 
target  of  thickness  1.5  cm. 
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These  initial  conditions  are  illustrated  in  the  preceding 
sketch.  Equation-of -state  and  constitutive  equations  are 
specified  in  (b)  and  (C)  below. 

2 .  Test  Problem  No.  2 

This  problem  differs  from  the  preceding  one  only  in 
regard  to  the  target,  which  becomes  a  1.0-cm-thick  aluminum 
plate,  as  sketched  below. 
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B.  MATERIAL  EQUATIONS  OF  STATE;  THE  HYDROSTATIC  STRESS 


The  equ:?tion  of  state  of  the  interacting  materials  is 
assumed  to  be  of  the  form 


P»(a  +Ay+Ay^+A  y’)(l  -  A  y) 

(1) 

+  (b  +  By  +  By*  +  B  y^^E 

where  the  A^  and  B^  are  constants,  P  is  pressure,  E  is 

specific  internal  energy  and  y  =  (p/Pj  '  ^)  is  a  measure  of 

compression  in  terms  of  the  density  p  and  initial  density 

p  .  This  form  is  sufficiently  general  for  representing  metals 
0 

over  the  pressure  ranges  of  interest. 

Values  of  the  material  constants  for  iron  and  aluminum 
are  listed  in  Table  1. 

C.  THE  DEVIATORIC  STRESS 


The  specified  relation  for  the  deviatoric  contribution 
to  the  stress  is  that  of  an  elastic-perfectly  plastic  material 
for  which  the  shear  yield  strength  K  satisfies  the  condition 


*=  ■  E  ‘ 

K  ■  0  for  E  ^  Ej^ 


(2) 


i.e.,  reduction  in  the  yield  stress  due  to  heating  is  accounted 
for  by  the  simple  linear  term  (l  - 

These  conditions  for  the  deviatoric  stress  are  easily 
formulated  along  lines  given  in  standard  texts  Starting  with 
the  usual  definition  of  the  strain  rate  deviator  tensor 
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TABLE  1 

Iron 

Aluminum 

A 

0 

0. 

0. 

A 

1 

1.236  X  10^2 

0.765  X  10^2 

A 

2 

2.452  X  10^2 

1.659  X  10^2 

A 

3 

5.139  X  io‘2 

0.428  X 

A 

4 

0. 

0. 

B 

0 

0. 

0. 

B 

1 

17.16 

5.94 

B 

2 

0  . 

0. 

B 

3 

0. 

0. 

C  (g/cc) 

0 

7.8 

2.79 

(mbar) 

0.006 

0.003 

(ergs/g) 

3.0  ^  10‘® 

7.0  X  10’ 

G  (mbar) 

1.93 

0.274 
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one  has  in  Cartesian  notation 


(3) 


For  the  elastic  regime  stress  deviator  increments  At 

'  J  w 

are  determined  from  the  relation  for  elastic  deformation 


2G  e, 


ij 


(4) 


where  G  is  the  rigidity  modulus.  This  elastic  relation  is 
used  to  determine  stress  increments  whenever  the  condition 


o 


ij 


a.  .  <  2K^ 


(S) 


is  satisfied.  When,  however,  the  stresses  given  by  (4)  would 
cause  (5)  to  be  violated,  one  instead  assumes  perfectly  plastic 
behavior  for  which  the  stresses  satisfy  the  von  Mises  condition 


(6) 


and  the  stress  increments  are  determined  from  the  Prandtl -Reuss 
equations 


o 


ij 


2G 


°k£  °kjl  _  \ 

-  0 .  .  I  . 

2K^ 


(7) 


It  is  easy  to  show  (multiplying  (7)  by  o—  and  performing  the 
indicated  summation)  that  the  use  of  (7)  to  compute  stress  in¬ 
crements  implies  that  o^j  Oj^j  *  0,  or  <^ij)  "  Hence 

(7)  is  compatible  with  (6)  in  the  sense  that  further  deformation 
of  the  plastic  type  causes  (6)  to  continue  to  be  satisfied. 
(Application  of  (4)  would  cause  (5)  to  be  violated.) 

Values  of  the  material  constants  G,  K  ,  Ej^  for  iron 
and  aluminum  are  given  in  Table  1. 
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The  above  elastic-perfectly  plastic  constitutive  equa¬ 
tions  are  idealized  versions  o£  the  material  behavior.  Such 
effects  as  the  dependence  of  stresses  on  strain  rates,  and 
the  increase  of  yield  due  to  work  hardening,  have  been  ignored 
in  these  test  problems  in  the  interest  of  simplicity. 

A  detailed  description  of  how  the  above  equations  are 
included  in  the  computer  program  will  be  given  when  the  HELP 
code  development  is  reported  as  part  of  Task  (1)  . 

D.  THE  FAILURE  CONDITION 


The  material  is  assumed  to  fail  only  due  to  hydrostatic 
tension.  Specifically,  the  iron  is  assumed  to  fail  at  a  rela¬ 
tive  volume  “  1*03  and  the  aluminum  fails  at 

(v/Vq)  »  1.015.  The  assumed  failure  thresholds  correspond  to 
tensions  of  approximately  35  kilobars  and  11  kilobars,  for 
unheated  iron  and  aluminum  respectively. 

More  sophisticated  failure  models  for  cavitation  or 
spalling  could  be  introduced.  During  the  past  few  years,  a 
considerable  amount  of  plane-wave  research  has  been  done  at 
Sandia  Corporation,  Los  Alamos,  General  Motors,  and  other 
laboratories,  to  determine  the  conditions  for  spallation 
failure  in  materials  of  interest.  A  summary  discussion  of  ex¬ 
perimental  work  and  proposed  mathematical  criteria  for  spalla¬ 
tion,  together  with  plane  wave  calculations  to  test  aspects  of 
these  models,  is  given  in  a  recent  S’  document.'-  ^  It  would 
not  be  difficult  to  incorporate  a  modern  time-dependent  spall 
criterion,  such  as  the  cumulative  damage  model^^^  by  Tuler 
and  Butcher,  into  a  two-dimensional  calculation.  Such  a  model 
would  cause  the  material  to  fail  at  a  tensile  level  which  de¬ 
pends  upon  the  local  tensile  history.  The  material  would 
thereafter  fail  to  support  tensions,  as  in  the  simple  instan¬ 
taneous  failure  model. 
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A  general  two-dimensional  treatment  of  the  material 
failure,  in  which  one  predicts  crack  propagation  directions 
and  velocities,  has  been  the  subject  of  some  attention  in 
the  present  program,  though  additional  work  would  be  required 
to  complete  and  test  such  a  capability. 

E.  REQUESTED  INFORMATION 

For  each  problem  plots  at  0.5  ysec  intervals  were 
requested  of  the  following  information: 

1.  Total  stress  (axial  direction)  versus 
Y  along  axis  of  symmetry. 

2.  Axial  velocity  versus  Y  along  axis  of 
symmetry . 

V  /  P  \ 

3.  Relative  volume  V“  \  ^ /  radial 

distance  at  five  levels  in  the  target: 

fa)  Target  front  surface 

(b)  One-quarter  way  through  target 

(c)  One-half  way  through  target 

(d)  Three-quarters  of  the  distance 
through  the  target 

(e)  Back  surface  of  the  target 

4.  Specific  internal  energy  at  these  positions. 

5.  Instantaneous  plastic  strain  rate  at  these 
positions . 

6.  Total  (integrated)  plastic  work  at  these 
positions . 

These  results  were  to  be  supplied  for  a  time  duration  of 
5  psec  in  Problem  No.  1,  and  for  3  ysec  in  Problem  No.  2. 
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In  regard  to  Items  (3)  through  (6),  the  requested 
quantity  (relative  volume,  specific  internal  energy,  etc.) 
is  plotted  versus  radial  distance  x,  along  that  locus  of 
moving  mass  elements  which  was  originally  (t*0)  at  the  re¬ 
quested  level  in  the  target.  Thus  the  plots  refer  to  the 
original  mass  elements  and  are  directly  comparable  to  Lagrangian 
plots  for  the  same  problem. 
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III.  NUMERICAL  SOLUTIONS  TO  THE  TEST  PROBLEM 
A.  GENERAL  REMARKS 


Solutions  have  been  obtained  to  20  psec  in  Test  Problem 
No.  1  and  to  10  ysec  in  Test  Problem  No.  2  in  addition  to  the  5 
and  3  psec  physical  times  which  were  requested.  Selected  re¬ 
sults  will  be  presented  for  the  extended  times,  but  the 
requested  detailed  plots  for  comparison  purposes  are  given 
only  to  5  and  3  ysec. 

No  special  difficulties  were  encountered  in  treating  the 
later  times  and  the  calculations  were  stopped  only  to  save  com¬ 
puting  costs. 

A  motion  picture  film  has  been  made  of  automatic  plots 
from  Test  Problem  No.  2,  showing  the  projectile-target  con¬ 
figurations  as  the  projectile  penetrates  the  target.  The  last 
frame  of  the  film  is  given  in  Section  III-C  as  Fig.  22. 

B.  RESOLUTION  AND  COMPUTING  TIME  INFORMATION 
1 .  Space  and  Time  Resolution 

For  both  of  the  problems  space  was  divided  into  cells 
of  dimensions  aZ  =  0.05  cm  and  aR  =  0.075  cm,  except  at 
large  radii  where  aR  was  increased  to  save  computing  costs 
in  remote  regions  of  the  target.  This  zoning  allowed  five 
cells  across  the  radius  of  the  projectile  in  either  problem, 
and  20  cells  across  the  target  in  Problem  No.  1. 

The  time  step  per  cycle  was  taken  to  be  the  minimum 
value,  for  any  cell  in  the  mesh,  of  the  quantity  0.4  AX/(c+u), 
where  AX  is  the  cell  dimension,  c  is  sound  speed  and  u 
is  the  absolute  magnitude  of  the  velocity  components  in  the 
AX  direction.  Thus  the  problems  were  run  at  0.4  times  the 
maximum  time  step  which  is  compatible  with  numerical  stability. 
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2 .  Computing  Time 

Detailed  computing  time  studies  were  made  of  Test 
Problem  No.  2.  At  late  times  in  the  interaction,  the  com¬ 
puting  time  per  average  (short  edit)  cycle  was  4.8  seconds. 
This  time  is  partitioned  among  the  major  subroutines  as 
follows . 


CDT 
CD  IT 
PHASE  1 
PHASE  3 
INFACE 

PHASE  2 


0.3  sec  (Equation  of  State  and  Time  Step") 

0.2  sec  (Energy  Conservation  Check  and  Edit) 

0.4  sec  (Pressure  Effects) 

1.0  sec  (Strength  Effects) 

1.5  sec  (Tracer  Motion,  plus  calculation 
of  mass  transport  from  mixed  and 
free  surface  cells) 

1.4  sec  (Transport) 


The  additional  time  to  treat  this  two-material  flow,  as  com¬ 
pared  to  the  analogous  one-material  problem  run  with  RPM,  is 
estimated  at  2/3  the  INFACE  time.  Thus  the  interface  and 
free  surface  treatments  increase  the  total  cost  of  the  cal¬ 
culation  by  about  20  percent.  This  percentage  increase  is  not 
fixed  from  one  problem  to  another,  however,  and  can  be  sub¬ 
stantially  different  for  problems  which  contain  a  substantially 
different  ratio  of  interface  cells  (free  surface  or  multi¬ 
material)  to  pure  cells. 

The  above  4.8  seconds  per  cycle  implies  8.5  minutes  for 
the  time  required  to  extend  the  problem  time  by  3  ysec.  The 
actual  computing  time  for  the  first  3  ysec  was  8.0  minutes, 
including  considerable  time  in  long  edits.  (This  time  is 
slightly  shorter  than  that  for  a  3  ysec  extension  at  later 
times,  because  the  entire  grid  is  not  yet  active  at  the 
earlier  times.) 
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Such  detailed  timing  studies  are  not  available  on 
Problem  No.  1  but  the  total  computing  time  for  the  first 
5.0  psec  was  11.0  minutes  (i.e.,  by  interpolation,  6.6  minutes 
for  3.0  psec).  This  'problem  is  somewhat  simpler  (no  material 
interface  nor  pressure  iteration  for  mixed  cells)  than  Problem 
No.  2  and  therefore  requires  less  computing  time. 

The  above  timing  studies  were  made  on  a  UNIVAr  1108 
computer. 

Economy  of  operation,  and  the  associated  freedom  to 
perform  parameter  studies  with  modest  computing  costs,  is 
one  of  the  strong  features  of  Eulerian  hydrodynamics. 

C.  PROJECTILE -TARGET  CONFIGURATIONS  VERSUS  TIME 
1 .  Test  Problem  No.  I 

Automatic  plots  of  the  pro jectile- target  configura¬ 
tions,  at  one-half  microsecond  intervals,  are  given  for  Test 
Problem  No.  1  as  Figs.  1  to  11.  Only  one-half  cf  the  inter¬ 
action  is  shown,  the  left  border  being  the  axis  of  symmetry. 

These  plots  are  constructed  using  the  positions  of  the 
moving  tracer  particles.  The  tracer  particle  interior  to 
the  target  are  plotted  as  points.  These  tracers  are  passive 
in  the  sense  that  they  do  not  "interact"  with  the  Eulerian 
solution  but  merely  move  with  the  local  fluid  velocity.  For 
the  projectile,  the  interior  tracers  are  connected  by  straight 
line  segments  and  are  also  passive  in  the  sense  described 
above.  The  interior  tracers  were  located  initially  in  the 
lower  left  corner  of  every  block  of  four  cells.  Their  loca¬ 
tions  at  any  time  gives  a  quick  impression  of  the  displace¬ 
ments  and  distortion  which  the  materials  have  undergone. 

Interface  tracers,  initially  five  per  Eulerian  cell 
dimension,  denote  the  position  of  the  moving  free  surfaces. 

(In  the  plots  the  free  surfaces  are  obtained  by  connecting 
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adjacent  tracers  with  straight  line  segments.)  These  inter¬ 
face  tracers,  which  are  also  advanced  with  the  velocity  of  the 
local  fluid,  are  not  passive  since  they  actively  affect  the 
basic  Eulerian  solution.  Transport  to  and  from  free  surface 
cells,  for  example,  takes  account  of  where  the  free  surface 
intersects  the  cell  boundaries;  further,  cells  are  evacuated 
when  the  surface  leaves  and  empty  cells  receive  no  mass  trans¬ 
port  until  entered  by  a  free  surface.  The  tracers  denoting 
the  projectile-target  interface,  in  this  like-material  impact, 
are  passive  since  no  particular  treatment  of  the  cells  con¬ 
taining  this  interface  is  required. 

Cells  which  have  spalled,  due  to  the  prescribed  failure 
condition,  are  marked  in  as  dark  rectangles.  Thus  peripheral 
rarefactions  have  caused  some  cavitation  at  the  interface  at 

1.5  microseconds.  And  the  projectile  opens  internally  near 
the  axis  between  2.0  and  3.0  usee,  only  to  close  again  before 

3.5  ysec.  In  this  latter  case,  the  oncoming  projectile  has 
filled  the  void  caused  by  the  rarefaction  from  the  sides  of 
the  projectile. 

2 .  Test  Problem  No.  2 

Similar  plots  of  the  iron-aluminum  impact  are  seen  as 
Figs.  12  through  22.  Here  the  plots  12-18  show  the  interaction 
at  one-half  microsecond  intervals  to  3.0  ysec,  and  times  5, 

7.5  and  10  ysec  are  shown  in  subsequent  plots.  The  plots 
12-21  show  the  fixed  Eulerian  grid  with  superimposed  contours 
of  the  moving  surfaces  and  the  projectile  target  interface. 
Figure  22  omits  the  grid  and  depicts  the  positions  of  the 
array  of  passive  tracers  which  were  originally  at  cell  centers 
within  the  projectile  and  target.  This  figure  shows  severe 
distortion  of  both  projectile  and  target  near  the  interface 
and  at  the  target  lip. 


14 


1 


3SR-201 


In  this  two-material  interaction  the  tracer  particles 
which  denote  the  projectile-target  interface  have  an  active 
role  in  the  Eulerian  solution.  In  the  calculation,  the  mass, 
momentum  and  energy  transports  (associated  with  the  moving 
steel  and  aluminum)  take  account  of  the  interface  position 
(obtained  by  interpolating  between  tracers) . 

Spalled  cells  are  darkened  in  the  plots.  It  is  seen 
that  spall  begins,  again,  near  the  interface  at  the  projec¬ 
tile  edge,  a  process  which  starts  at  1.0  ysec  and  increases 
throughout.  (The  projectile  does  not  fail  near  the  axis,  as 
in  Problem  No.  1.)  Lack  surface  target  spall  is  evident, 
beginning  at  2.0  psec,  reaching  a  maximum  at  2.5  to  3.0  ysec, 
and  diminishing  somewhat  at  later  times  as  the  target  material 
begins  to  overtake  the  spalled  material.  This  back  surface 
spall  is  also  evident  (as  tracer  particle  separation)  in  the 
motion  picture  which  was  referred  to  in  (A)  above. 

D.  PLOTS  OF  REQUESTED  INFORMATION 
1.  Problem  No.  1 

(a)  Total  Stress  (Axial  Direction) 

Versus  V  Along  Axis  of  Symmetry 

These  results  are  plotted,  for  half-microsecond  inter¬ 
vals  to  5  ysec,  as  Figs.  23  through  32.  In  all  plots  the 
initial  position  of  the  interface  is  at  1.4  cm  and  distance 
into  the  target  increases  to  the  right.  At  the  earliest  time, 
0.5  ysec,  nearly  symmetrical  150  kbar  shocks  have  propagated 
into  both  the  target  and  projectile.  At  1.0  psec  these  shocks 
have  attenuated  to  approximately  half  this  strength  due  to 
shock  divergence  and  to  rarefactions  from  the  free  surfaces. 
By  1.5  ysec,  a  substantial  tensile  region  (20  kbars)  is  seen 
on  the  axis;  this  tensile  region  is  predominantly  in  the  pro¬ 
jectile  material,  near  the  projectile- target  interface.  The 
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plot  at  2.0  psec  shows  the  effect  of  cavitation  failure  on 
the  projectile  axis,  occurring  between  1.5  and  2.0  usee  when 
tensions  exceed  some  35  kbars  (see  Section  II-2  for  the 
assumed  failure  condition).  This  tension  results  in  a  rare¬ 
faction  into  the  neighboring  material,  and  slight  noise  in 
the  solution  due  to  the  sudden  increase  of  pressure  to  zero 
when  the  failure  condition  is  satisfied.  Tensions  are  also 
evident  in  the  rear  of  the  projectile  at  2.5  psec,  but  these 
are  seen  later  to  be  insufficient  to  cause  cavitation.  By 
3.5  psec,  further,  the  oncoming  projectile  has  caused  the 
previously  spalled  cells  to  again  be  compressed,  and  positive 
pressures  are  seen  at  all  positions  on  the  axis.  At  5  psec 
the  stresses  are  compressive  and  the  maximum  stress,  approxi¬ 
mately  20  kbars,  is  in  a  region  approximately  one-fourth 
through  the  target. 

(b)  Velocity  Versus  Y  Along  Axis  of  Symmetry 

The  axial  velocity  plots  are  seen  as  Figs.  33  through 
43.  These  plots  are  essentially  self-explanatory.  The  results 
are  believed  to  be  realistic,  except  for  the  expected  smearing 
of  the  shock  front  over  several  cells  and  the  observed  small 
oscillations  in  the  velocity  profile.  It  can  be  observed 
that,  at  the  latest  time  of  5  psec,  the  back  of  the  projectile 
has  decelerated  only  about  six  percent  in  velocity,  from  the 
original  9  x  10**  cm/sec  to  8.4  x  lo**  cm/sec.  (At  this  time, 
the  energy  in  the  problem  remains  57  percent  kinetic,  as  indi¬ 
cated  by  the  energy  values  given  at  the  bottom  of  Fig.  11.) 

(c)  Relative  Volume  Versus  Radial  Distance  at 
Various  Levels  in  the  Target 

Plots  of  relative  volume  X-  (■  are  seen  as 

V  j  '  Y  P ' 

Figs.  44  through  53.  Compressions  to  ^  *  .9  are  seen  at 
the  earliest  time  (0.5  psec)  at  the  target  front  surface  in 
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the  vicinity  of  the  axis.  By  1.0  ysec  the  compression  in 
this  region  has  diminished  by  about  a  factor  of  five  ®  *98^ 
and  by  1.5  ysec  the  compression  has  become  a  distension 
(■^  -  1.02),  as  seen  in  Fig.  46.  This  distension  corresponds 
to  the  tension  noted  previously  (Part  (a)  above)  in  connection 
with  the  results  for  stress  as  a  function  of  distance  along  the 
^is  of  symmetry.  By  times  of  2.5  to  5  ysec,  small  oscillations 
are  evident  in  these  plots  of  relative  volume  versus  radial 
distance,  especially  at  the  target  front  surface.  These 
oscillations  are  due  primarily  to  the  sudden  discontinuity 
which  occurs  in  the  interaction  whenever  the  material  in  a 
given  cell  reaches  its  failure  threshold.  (Somewhat  similar 
effects  were  seen  previously  in  the  axial  pressure  plots.) 

To  good  accuracy  at  the  stress  levels  in  this  problem,  these 
plots  of  relative  volume  are  equivalent  to  pressure  plots  at 
the  same  times  and  levels,  since  pressure  is  proportional  to 
compression  for  small  compressions  of  the  cold  material. 

(d)  Specific  Internal  Energy  Versus  Radial 
Distance  at  Various  Target  Levels 

These  results  are  given  as  Figs.  54-63.  The  greatest 
specific  internal  energies,  around  1.25  x  10*  ergs/g,  occur 
at  early  times  (see  Fig.  54,  the  0.5  ysec  plot)  and  are  simply 
the  specific  internal  energy  of  the  initially  shocked  material. 
This  particular  energy  (that  due  to  the  plane  shock  at  early 
times)  is  very  easily  determined  theoretically  from  the  shock 
relation  AE  ■  JjUp  where  the  shock  particle  velocity  Up, 
from  symmetry,  must  be  one-half  the  9.0  x  10®  velocity  of  impact 
The  resulting  energy  is  1.01  x  10®  ergs/g.  The  agreement  is 
reasonable,  the  difference  between  this  exact  value  and  those 
in  Fig.  54  being  due  to  transient  oscillations  in  the  early 
stages  of  the  numerical  solution. 
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At  the  later  times,  the  energies  are  seen  to  diminish 
as  pressures  are  relieved,  and  the  material  expands  from  its 
shocked  state.  Once  the  pressures  have  attenuated  to  low 
stress  levels,  as  at  the  5  psec  time  depicted  in  Fig.  63, 
the  observed  specific  internal  energies  are  attributable  to 
a  combination  of  (a)  shock  heating,  and  (b)  work  done  against 
the  plastic  forces  (the  additional  mechanism  of  energy  change, 
adiabatic  expansion,  being  negligible  at  low  stress  levels.) 

As  seen  in  the  late  time  plots,  the  specific  internal 
energies,  near  the  target  front  surface,  peak  at  a  radius  of 
.3  or  .4  cm,  so  that  there  is  a  (local)  minimum  of  specific 
internal  energy  at  the  axis.  This  apparent  anomaly,  since 
one  might  expect  a  maximum  at  the  axis  of  symmetry,  has  been 
the  subject  of  a  simple  analysis.  The  effect  arises  because 
the  work  done  against  the  plastic  forces  (plotted  subsequently 
as  Figs.  74  through  83)  is  a  maximum  for  the  high-shear  region 
near  the  projectile  edge  (at  radius  .375  cm)  and  therefore  has 
a  relative  minimum  at  the  axis.  Further,  a  quantitative  check 
can  be  made  on  the  value  of  the  specific  internal  energy  at  the 
center  of  impact  (target  front  surface  at  the  axis  of  symmetry) 
the  shock  heating  is  given  by  the  energy  on  the  Hugoniot 
(calculated  as  1.01  x  lo^  ergs/g  in  the  preceding  paragraph) 
minus  the  loss  in  energy  /pdV  due  to  the  adiabatic  expansion. 
The  latter  integral  has  been  evaluated  using  a  special  equation 
of-state  program,  giving  0.88  x  lo*  ergs/g.  Thus  the  net  shock 
heating  is  1.3  x  lo*  ergs/g.  The  work  per  unit  volume  against 
the  plastic  forces  (Fig.  83)  is  seen  to  be  W  ■  4.5  x  io‘* 
megabars  «  4.5  x  lo’  ergs/cm*;  the  work  per  gram  is  therefore 
Ep  ■  W/p  *  (4.5  X  10®/7.85)  ergs/g  ■  5.7  x  10*  ergs/g.  The 
sum  of  the  irreversible  shock  heating  and  this  plastic  work 
is  therefore  E  =  (1.3  +  5.7)  10*  ergs/g  ■  7.0  x  10*  ergs/g. 
This  value  of  the  specific  internal  energy  at  late  times  is 
in  very  good  agreement  with  the  result  7.5  x  10*  ergs/g  seen 
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in  Fig.  63,  and  which  was  obtained  by  integrating  the  hydro¬ 
dynamic  equations  of  motion.  From  this  discussion  we  see 
then  that  (a)  the  commuted  value  of  the  specific  internal 
energy  at  the  impact  center  is  reasonable,  and  (b)  the  in¬ 
crease  of  specific  internal  energy  with  radial  distance,  along 
the  target  front  surface,  is  attributable  to  the  large  amount 
of  work  against  plastic  forces  in  the  region  of  maximum  target 
shear. 


(e)  Instantaneous  Plastic  Strain  Rate  Versus 
Radial  Distance  at  Five  target  leveTs 

These  results  are  seen  as  Figs.  64  through  73. 

The  plastic  strain  rates  data  represent  the  least 
smooth  of  the  results  to  be  presented  because  these  data  are 
intensive  quantities  obtained  by  differentiation  of  other 
intensive  quantities  (the  velocities)  with  respect  to  the 
space  variables.  Small  oscillations,  e.g.,  those  at  the  1/4 
level  at  1.0  ysec,  are  apparently  spurious,  but  the  data  are 
reliable  in  regard  to  general  trend  and  the  time -averaged 
(or  integrated)  results. 

Very  high  strain  rates  are  evident  at  early  times  in 
the  vicinity  (r  ■  0.375  cm)  of  the  projectile  edge,  corres¬ 
ponding  to  the  maximum  shear  region.  By  the  later  times  this 
region  of  maximum  shear  occurs  at  larger  radii  (e.g.,  at 
r  *  6  cm  in  the  5  ysec  plot.  Fig.  73).  Of  the  various  levels 
in  the  target  the  highest  strain  rates  are  always  for  the 
material  which  was  originally  near  the  front  target  surface. 

(f)  Plastic  Work  Versus  Radial  Distance 
at  Various  Levels  in  the  TargFt 

Plots  of  the  total  plastic  work  per  unit  volume  of 
material  are  presented  as  Figs.  74-83.  The  results  are  essen¬ 
tially  self-explanatory  and  show,  as  was  indicated  by  the 
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strain-rate  results  described  in  (d)  above,  that  most  of  the 
plastic  work  is  concentrated  in  the  impact  region  near  the 
target  front  surface.  Further,  as  was  pointed  out  above,  the 
plastic  work  is  somewhat  higher  at  the  projectile  radius 
distance  than  at  r  =  0,  causing  a  relative  minimum  in  the 
specific  internal  energy  at  the  axis.  At  5.0  ysec  the  inte¬ 
grated  values  of  the  total  plastic  work  for  the  entire 
problem  (from  a  separate  calculation)  is  4.2  x  lO’  ergs.  This 
compares  with  a  total  internal  energy  of  6.8  x  10*  ergs.  Thus, 
work  against  the  plastic  forces  has  accounted  for  some  62  per¬ 
cent  of  the  gain  in  internal  energy. 

2 .  Problem  No.  2 

(a)  Total  Stress  (Axial  Direction) 

Versus  V  Along  Axis  of  Symmetry 

These  plots  of  total  stress  are  seen  as  Figs.  84  through 
89.  At  the  earliest  time,  0.5  psec,  125  kbar  shocks  are  seen 
to  propagate  in  both  directions  from  the  impact  surface.  As 
a  result  of  rarefactions  and  shock  divergence,  these  shocks 
attenuate  to  75  kbar  by  1.0  ysec.  Subsequent  plots  show  con¬ 
tinuing  shock  attenuation,  accompanied  by  oscillations  near 
the  interface  (Fig.  86,  time  1.5  ysec)  and  near  the  back  sur¬ 
face  of  the  target  (Fig.  87,  time  2.0  ysec).  These  fluctua¬ 
tions  in  the  stresses  correlate  with  material  failure  in 
tension  (see  Figs.  12-22)  which  begins  at  the  interface  and 
back  surface  just  before  the  times  of  1.5  and  2.0  ysec 
respectively.  This  sudden  release  of  tension,  occurring 
when  the  material  satisfies  the  prescribed  failure  conditions 
(see  Section  II-D)  gives  rise  to  local  disturbances. 

(b)  Velocity  Versus  Y  Along  Axis  of  Symmetry 

These  plots  for  Problem  No.  2  are  seen  as  Figs.  90 
through  95.  Two  effects  will  be  noted.  At  2.0  ysec.  Fig.  93, 
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the  projectile  material  in  the  vicinity  of  the  interface 
actually  has  a  velocity  which  is  slightly  greater  than  the 
initial  velocity  of  impact,  by  a  few  percent.  Whether  this 
apparent  slight  acceleration  is  real  is  not  quite  clear. 

Also,  at  the  latest  times,  spallation  is  evident  near 
the  back  surface  of  the  target,  where  the  material  near  the 
back  surface  is  seen  to  have  a  greater  velocity  than  the  tar¬ 
get  material  nearer  to  the  projectile.  This  evidence  of 
spallation  correlates  well  with  spallation  results  presented 
in  Figs.  12-22. 


Relative  Volume  Versus  Radial  Distance 
at  Various  Levels  in  the  Target 


These  plots,  for  the  six  requested  times,  are  shown  as 
Figs.  96  through  101.  At  the  earliest  time  (.5  ysec,  Fig.  96) 
the  maximum  compressions  ^  seen  to  occur 

near  the  target  front  face,®i.e.,  at  the  front  surface  level 
and  the  one-fourth  level.  At  1.0  ysec  (Fig.  97)  maximum 
compression  occurs  in  the  middle  of  the  target,  as  seen  by 
comparing  the  third  plot  with  the  other  four.  At  still  later r 
times  (e.g.,  1.5  ysec,  Fig.  98)  the  compression  pulse  has  . 
advanced  further,  so  that  maximum  compression  occurs  at  the 
three-quarter  station.  The  behavior  is  indicative  of  a  pressure 
pulse  which  propagates  at  the  aluminum  shock  wave  velocity. 


The  occasional  irregularities  in  the  plots  at  the 
target  front  surface  level  (vicinity  of  fourth  or  fifth  point 
in  Fig.  97,  vicinity  of  fifth  or  sixth  point  in  Fig.  101) 
reflect  the  fact  that  these  points  are  in  the  (small)  region 
of  the  highly  deformed  crater  lip,  where  the  material 
characteristics  are  rapidly  changing. 
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(d)  Specific  Internal  Energy  Versus 

Radial  Distance  at  Various  Target  Levels 

These  results,  for  the  various  times,  are  seen  as 
Figs.  102  through  107  inclusive.  As  might  be  expected,  specific 
internal  energies  are  highest  near  the  target  front  surface 
and,  also,  tend  to  be  somewhat  higher  at  the  projectile  radius 
distance  than  at  the  axis  or  at  greater  radii. 

(e)  Instantaneous  Plastic  Strain  Rate  Versus 
Radial  Distance  at  Five  Target  Levels 

These  data  for  Test  Problem  No.  2  are  displayed  as 
Figs.  108  through  113  inclusive. 

(f )  Plastic  Work  Versus  Radial  Distance  at 
Various  Levels  in  the  Target 

Plots  of  the  total  plastic  work  are  shown  in  Figs.  114 
through  119.  The  earliest  plots  exhibit  a  concentration  of 
plastic  work  in  the  high  shear  region,  near  the  target  front 
face  at  the  projectile  radius.  Eventually,  this  effect  is 
over-ridden  (see  3.0  ysec.  Fig.  119)  as  much  more  plastic  work 
is  performed  within  the  target.  In  this  impact  the  softer  and 
lighter  aluminum  appears  to  be  flowing  around  the  iron  projectile, 
causing  an  associated  build-up  within  the  aluminum  of  work  done 
against  plastic  forces.  This  process  continues  long  past 
3  (jsec,  at  which  time  separate  edits  of  kinetic,  internal 
and  plastic-work  energies  show  that  the  energy  in  the  problem 
is  still  82%  kinetic,  with  only  11%  of  the  energy  having 
been  expended  against  plastic  forces.  At  the  later  time  of 
10.0  psec  (Fig.  22  shows  the  configuration  at  this  time)  the 
energy  is  53%  kinetic,  with  nearly  all  of  the  remaining  energy 
being  lost  due  to  work  against  plastic  forces. 
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IV.  CONCLUDING  REMARKS 

The  preceding  pages  consist  of  the  information  requested 
by  the  JTCG/ME.  The  results  given  are  complete  and  have  been 
presented  in  such  a  form  (at  locations  of  moving  tracer  points) 
that  comparisons  can  be  made  directly  with  the  corresponding 
Lagrangian  code  results . 

Much  related  work  has  been  described  in  previous  reports 
This  includes  the  SLUG  and  DORF  code  calculations  mentioned  in 
the  Introduction,  and  those  times  in  the  HELP  calculations  that 
are  later  than  those  requested.  Also,  a  number  of  exploratory 
calculations  were  made  to  evaluate  the  effects  of  smaller  time 
steps  or  slight  modifications  in  the  transport  from  mixed  cells 
These  latter  calculations  did  not  differ  significantly  from 
those  presented  here. 

Much  additional  information  pertaining  to  the  above 
calculations  can  be  obtained  from  the  printed  numerical 
solutions.  S*  will  be  happy  to  supply  such  information  in 
response  to  queries  or  can  loan  the  printed  results  to  respon¬ 
sible  parties. 

In  regard  to  the  ongoing  program,  it  is  expected  that 
Item  (2)  of  the  Introduction,  the  application  of  Eulerian 
methods  to  give  the  most  realistic  possible  descriptions  of 
armor  penetration  processes,  will  receive  primary  attention 
in  the  months  immediately  ahead. 
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TEST  PROBLEM  NO.  1 

Projectile-Target  Configuration  with  Spalled  Cells  Shaded 

For  T  »  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0, 

4.5,  and  5.0  ysec. 

(In  target  interior  tracers  plotted  as  points) 

(In  projectile  interior  tracers  connected  by  straight 
line  segments) 

(Interior  tracers  originally  located  at  bottom  left  corner 
of  every  block  of  four  cells) 
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TEST  PROBLEM  NO.  2 


Proiectile-Target  Configuration  Superimp 
on  Eulerian  Grid  with  Spalled  Cells  Shad 
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For  T 


0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  5.0, 
7.5,  and  10.0  ysec. 
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test  problem  no.  1 


Total  Stress  fAxial  pirectionl  Versus — Y_ 
Along  Axis  oi:  Symmetry 

For  T  -  0.0,  0.5,  1.0,  1.5,  2.0.  2.5,  3.0, 
4.5,  and  5.0  ysec. 


3.5,  4.0, 


PB  ■  Projectile  Back 
PTI  ■  Projectile  Target  Interface 
TB  "  Target  Back 
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TEST  PROBLEM  NO.  1 

Axial  Velocity  Versus  Y  Along  Axis  of  Symmetry 

For  T  -  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0, 
4.5,  and  5.0  usee. 

PB  “  Projectile  Back 
PTI  =  Projectile  Target  Interface 
TB  ■  Target  Back 
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Figure  40 
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IRON  IMPACT  TEST  PRQB 

Figure  43 


3SR-201 


TEST  PROBLEM  NO.  1 

Relative  Volume  Versus  Radial  Distance  at 
Five  Locations  in  the  Target 

For  T  ■  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0, 
4.5,  and  5.0  ysec . 

A  *  Target  Front  Surface 
B  =  One-Quarter  Distance  into  Target 
C  =  One-Half  Distance  into  Target 
D  “  Three-Quarters  Distance  into  Target 
E  “  Back  Surface 
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X  CCM5  TIME-  1.0000-06 

Figure  45. 
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IRON  inPACT 


RH0Z/nn3  vs  X  (CMJ  TIMES  4.5000-06 


3SR-201 


TEST  PROBLEM  NO.  1 


For  1  ■  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0 
4.5,  and  5.0  ysec. 

A  ■  Target  Front  Surface 
B  ■  One-Quarter  Distance  into  Target 
C  *  One  Half  Distance  into  Target 
D  ■  Three-Quarters  Distance  into  Target 
E  ■  Back  Surface 


Specific  Internal  Energy  Versus  Radial  Distance 
at  Five  Locations  in  tne  Target 


IRON  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/G}  VS  X  CCM)  TIME-  S. 0000-07 

Figure  54. 


IR9N  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/G3  VS  X  CCM3  TIME«  1.0000-06 


IRCSN  -  IRON  IMPACT 


SPECIFIC  INTZ\.;.'.L  E.'iHF.SY  (Ei^GS/C)  VS  X  (CM) 

Figure  5b 
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SPECIFIC  INruroCAL  ENZRUY  CERGS/G)  VS  X  CCM]  TlriE=  4.0000-06 


IRON  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/G)  VS  X  CCM)  TIHE=  4.S0C0-06 

Figure  62. 


IRON  -  IRON  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/G3  VS  X  (CM)  TIME=  5.0000-06 
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TEST  PROBLEM  NO.  1 

Instantaneous  Plastic  Strain  Rate  Versus  Radial 
Distance  at  Five  Locations  in  the  Target 

For  T  ■  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0, 
4.5,  and  4.0  ysec. 

A  -  Target  Front  Surface 
B  -  One-Quarter  Distance  into  Target 
C  ■  One-Half  Distance  into  Target 
D  ■  Three-Quarters  Distance  into  Target 
E  ■  Back  Surface 
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PLASTIC  STRAIN  RATE  C1/SEC3  VS  X  tCM)  TIME=  1.5000-06 
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4.0000*0$ 


100 


IRON  IMPACT 


PLASTIC  STRAIN  RATE  Cl/SEC)  VS  X  CCM)  TIME=  A. 0000-06 


IRON  INPACT 
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PLASTIC  strain  RATE  tl/SEC)  VS  X  (CN)  TiriE=  A.SOOO-06 


3SR-201 


TEST  PROBLEM  NO.  1 

Plastic  Work  Versus  Radial  Distance  at  Five 
Locations  in  the  Target 

For  T  -  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  4.0, 
4.5,  and  5.0  ysec. 

A  ■  Target  Front  Surface 
B  -  One-Quarter  Distance  into  Target 
C  ■  One-Half  Distance  into  Target 
D  *  Three-Quarters  Distance  into  Target 
E  ■  Back  Surface 
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IRON  INPACT 


PLASTIC  WORK  (MBJ  VS  X  tCM)  TIME=  1.5000- 

Figure  76. 
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IRON  .  IRON  IMPACT 
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PLASTIC  WORK  tMBJ  VS  X  CCM)  TIt«s  3.5000-06 

Figure  80. 


IRON  IMPACT 
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TEST  PROBLEM  NO.  2 

Total  Stress  (Axial  Direction)  Versus  Y 
Along  Axis  of  Symmetry 

For  T  -  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  ysec. 

PB  “  Projectile  Back 
PTI  =  Projectile  Target  Interface 
TB  =  Target  Back 
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TtHE«  5. 0000-07  PROBs  2.6S000«-0I  CYCLE  NO 


3SR-201 


117 


ALUniNUn  IMPACT 
Figure  84  . 
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Figure  87 


time*  2.5000-06  PROS*  2.65000*^01  CYCLE  NO 
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IRON  -  ALUMINUM  IMPACT  -  HELP-5 


TIME*  3.0000-06  PROB*  2.6S000«^01  CYCLE  NO  107 
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IRON  -  ALUMINUM  IMPACT  -  HELP-5 


3SR-201 


TEST  PROBLEM  NO.  2 


Ax:al  Velocity  Versus 


Along  Axis  of  Symmetr) 


For  T  =  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  ysec 


PB  =  Projectile  Back 
PTI  =  Projectile  Target  Interface 
TB  =  Target  Back 
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TEST  PROBLEM  NO.  2 


Specific  Internal  Energy  Versus  Radial 
at  Five  Locations  in  tne  TargeT 

For  T  -  0.0,  0.5,  1.0,  1.5,  2.0,  2.5, 


Distance 


and  3.0  y 


A  =  Target  Front  Surface 
B  =  One-Quarter  Distance  into  Target 
C  =  One-Half  Distance  into  Target 


IRON  -  ALUHINUN  IHPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/GI  VS  X  CCM)  TIME=  5.0000-07 


IRON  -  ALUMINUM  IMPACT 


I 


SPECIFIC  INTERNAL  ENERGY  CERGS/GJ  VS  X  CCM)  TIME=  I.SOOO-06 

Figure  104. 


IRON  -  ALUMINUM  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CERGS/G)  VS  X  f^MJ  TIME=  2.0000-06 


ALUMINUM  IMPACT 


SPECIFIC  INTERNAL  ENERGY  CER6S/6)  VS  X  CCMJ  TIME=  2.S000-06 


ALUMINUM  IMPACT 


0000*0 
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TEST  PROBLEM  NO.  2 

Instantaneous  Plastic  Strain  Rate  Versus  Radial 
Distance  at  Five  Locations  in  the  TargeT 

For  T  =  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  ysec. 

A  =  Target  Front  Surface 
B  =  One-Quarter  Distance  into  Target 
C  =  One-Half  Distance  into  Target 
D  =  Three-Quarters  Distance  into  Target 
E  =  Back  Surface 
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ALUMINUM  IMPACT 


O 


PLASTIC  STRAIN  RATE  C1/SEC3  VS  X  CCMJ  TIME*  5.0000-07 

Figure  108. 
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PLASTIC  STRAIN  RATE  tl/SEC)  VS  X  tCM)  TIME*  2.5000-06 


PLASTIC  STRAIN  RATE  C 1/SEC)  VS  <  iCN)  TIHE*  3.0000-06 

Figure  113. 
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TEST  PROBLEM  NO.  2 

Plastic  Work  Versus  Radial  Distance  at  Five 
Locations  in  the  Target 

For  T  =  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  ysec. 

A  =  Target  Front  Surface 
B  =  One-Quarter  Distance  into  Target 
C  =  One-Half  Distance  into  'iarget 
D  =  Three-Quarters  Distance  into  Target 
E  =  Back  Surface 
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PLASTIC  WORK  (MB)  VS  X  (CM)  TiME=  2.5000-08 

Figure  118. 


